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A theory is developed for magnetically confined Fermi gas at low tem- 
perature based on the density functional theory. The theory is illustrated 
by numerical calculation of density distributions of Fermi atoms 40 K with pa- 
rameters according to DeMarco and Jin's experiment [Science, 285(1999)1703]. 
Our results are in good agreement with the experiment. To check the the- 
ory, we also performed calculations using our theory at high temperature and 
compared very well to the result of classical limit. 
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In a recent publication, DeMarco and Jin reported their work on Fermi atoms at low 
temperature in confined geometries. They employed an evaporative cooling strategy to 
cool a magnetically trapped gas of 7 x 10 5 40 K atoms to 0.5 of the Fermi temperature 
Tp. An Ioffe-Pritchard-type magnetic trap |2| provides a cylindrically symmetric, harmonic 
potential with an axial frequency of u z = In x 19.5 Hz and a variable radial frequency. 
The radial frequency, as well as the minimum magnetic field, can be smoothly varied from 
uj r = 2tt x 44 Hz to 2it x 370 Hz by changing the current in a pair of Helmholtz bias coils. 
In the experiment presented by DeMarco and Jin, a single-component gas of the F = 9/2, 
Trip = 9/2 atoms is produced, where F is the total atomic spin and mp is its magnetic 
quantum number. The time-of-flight image was taken by suddenly switching off the current 
that provides the magnetic trapping field, which allowed the gas to expand freely for 15 to 20 
ms. The absorption shadow, generated by illumination of the expanded gas, was imaged onto 
a charged-coupled device array. They detected the emergency of quantum degeneracy in a 
trapped gas of Fermionic atoms and observed a non-classical momentum distribution and 
found that the total energy of the gas is larger than the classical expectation. They observed 
not only the momentum distribution but also the confined gas itself, which approaches a 
fixed size as T approaches zero. 

In this letter, we present a general theory of the low-temperature density profile of an ideal 
Fermi gas trapped within an arbitrary potential well. Then we give the results of numerical 
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calculation to compare to the experiment. The approach taken here is based on the density 
functional theory (DFT) which was first introduced by Kohn and Hohenberg |J in 

the context of ground state energy of quantum systems. It was developed and extended to 
excited states and to finite temperatures by some researchers and has become the effective 
first principle calculational method for the electronic and structural properties of a large 
variety of condensed matter systems. Considering a single- component system within an 
external potential at T / OK, the DFT asserts that the true density distribution of a 
system in an external field V(r) is the one that leads to the minimum of the following 
functional, 

F[n(r)] = F [n(r)] + J [V(r) - //] • n(r)dr, (1) 

and the minimum is the the Helmholtz free energy of the system. Here /i is the chemical 
potential in the above expression. While Fo[n(r)] is a temperature dependent functional 
of n(r) only, the dependence of i<o[n(r)] on external potential V(r) is only through the 
dependence of n(r) on V(r). This is a very strong statement which means that if we know 
the functional form of F [n(r)] we can get all the physical quantities of interest. For a given 
external potential, the functional derivative of F[n(r)] with respect to n(r) should equals to 
zero: 
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This gives an equation for n(r), and when n(r) is obtained from the solution of Eq. (0), we 
substitute it back to (|1|) to get the free energy and other quantities can be obtained simply 
by differential. It is hard to get the functional form of -F [n(r)], various approximations are 
employed in practical calculations. One of the commonly used approximation is the local 
density approximation(LDA) which usually gives excellent results to real systems. We use 
here the LDA in our theory of confined Fermi systems. In this approximation the functional 
F [n(r)] is assumed to be 

F„[n(r)] = J dr/„(n(r)). (3) 

And we assume that the free energy density fo(n(r)) has also the same functional form as 
the free energy of an uniform system of density n, i.e, fo(n) is the Helmholtz free energy 
density of an uniform system with density n. Then we can obtain 

5n(r) dn{v) + [V [ } ^ [ ) 

According to Eqs. (fj) and (|J), we have 
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df (n)/dn =n~ V(r). 



(5) 



For the ideal Fermi gas, we know 



n = ^3/3/2(2), (6) 



and 



kT 

fo = —r^-f 5/2(2) + nkT\ogz } (7) 



here 
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where k is Boltzmann's constant, A = J 2tiTi 2 j mkT the thermal wave length, z = e^ the 
fugacity, and /xo the chemical potential of the uniform system. So we get 

df (n)/dn = kT log z. (8) 

The value of /x is determined from the normalization of n(r, /x): 

N = J n(r,n)dr (9) 

where we have indicated the /x dependence of n which follows from Eq. (Q). Once having 
determined /x from Eq. (||) we substitute it back into Eq. @ to get n(r). Oliva has developed 
an approximate form of the Helmholtz free energy of the uniform system for different regimes 
of density 0. However, in our theory such an expression of fo with respect to n is not 
required. 

Now we use the confining potential well, according to the experiment by DeMarco and 
Jin, to calculate the density distributions of Fermi atoms. 
We consider a harmonic potential form 

t r / \ 22 22 

V{r,z) = -muj r r + -muj z z 

which is cylindrically symmetric, where the axial frequency uj z = 2tt x 19.5Hz and the radial 
frequency u r = 2tt x 137Hz. The potential is independent of the azimuth angle. We have 
used the parameters given by the experiment setup, which is iV = 7 x 10 5 and T = Tp/2, 
with T F = U^ujzUJ^N) 1 ^ jk |8|||. The value of T F is 0.6/xK for a million atoms in the 
uj r = 2-7T x 137Hz trap. 
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FIG. 1. Atomic density as function of z with different values of r for 40 K atoms under mag- 
netic confinement. Total number of atoms N = 7 x 10 5 , (a) T = 0.5Tp, (b) T = 0.05Tp, with 
T F = 0.6fiK. 
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FIG. 2. Atomic density as function of r with different values of z for 40 K atoms under magnetic 
confinement. Other parameters are the same as Fig. 1. 

The particle density distributions are calculated by using the present method. Figures 
1 and 2 are plots of the density as function of z and r for two different temperatures. As 
observed in the experiment that the confined gas approaches a fixed size, we calculated the 
root-mean-square radius of the system and get Trms = 0.7Rf at T = 0.5Tp and trms = 
0.5-Rp at low temperature T = 0.05Tp, where Rp = J2kTp/rnoj^ [§,0- These results are 
coincident with the experiment, where they obtained trms — 0.6Rp as T approaches zero 
and r RMS = 0.9Rp at T/T F = 0.5 QTJ • The difference between theory and experiment may 
due to the systematic uncertainty in T/Tp within the experiment. To check the validity of 
our method we calculated the case for T/Tp = 2 with other conditions unchanged. At high 
temperature (i.e., in the classical limit T » Tp), we have ||: 
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/i=-A;Tlog[6(-^) 3 ] (10) 

F 

In our calculation the chemical potential of the system is /i — —5.86 x 10~ 29 J at T/T F = 2, 
while Eq. (10) gives —5.92 x 10~ 29 J, which is satisfactory. The comparison of density 
distributions of our theory with the classical statistics can be seen in Fig. 3. 
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FIG. 3. Atomic density as function of r with different values of z for 40 K atoms under magnetic 
confinement. Parameters are the same as Fig. 1 except T = 2Tp. The solid lines are calculations 
from DFT as described in this work and dashed lines are from the Boltzmann distributions. 

In summary we calculated the physical properties of the Fermi gas confined within a 
potential well at none-zero temperature. This method is based on the density functional 
theory under local density approximation. The case of 40 K atoms in a magnetic confinement 
is calculated and the result agrees with recent experiment. Here the interaction between 
atoms is completely neglected and the study of effect of interaction on the confined Fermi 
atomic gas is underway. Further studies of quantum degenerate Fermi gases in different 
confinement and interactions are also the important and interesting subjects. 
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